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Abstract 

The internal shocks scenario in relativistic jets is used to explain the variability 
of the blazar emission. Recent studies have shown that the magnetic field signif¬ 
icantly alters the shell collision dynamics, producing a variety of spectral energy 
distributions and light-curves patterns. However, the role played by magnetization 
in such emission processes is still not entirely understood. In this work we numer¬ 
ically solve the magnetohydodynamic evolution of the magnetized shells collision, 
and determine the influence of the magnetization on the observed radiation. Our 
procedure consists in systematically varying the shell Lorentz factor, relative ve¬ 
locity, and viewing angle. The calculations needed to produce the whole broadband 
spectral energy distributions and light-curves are computationally expensive, and 
are achieved using a high-performance parallel code. 


1 Introduction 

Blazars are a class of radio-loud active galactic nuclei (AGNs) whose jets are pointing 
near the line of sight to the observer Ifl3l . and are known for showing the most rapid 
variability of all AGNs. Extensive observations indicate that the jets of blazars are rel¬ 
ativistic. Their most remarkable characteristic is the presence of flares in the X- and 
y-ray bands, usually wih a duration of the order of a few hours. The internal shock (IS) 
model ESI, invoked to explain this variability 001321 , is an idealized model where an 
intermittently working central engine ejects shells of magnetized plasma which collide 
due to their velocity differences. As a consequence internal shocks are formed, ac¬ 
celerating particles at the shock fronts and producing the non-thermal, highly variable 
radiation. 

Our long-term objective for is the study of the influence of the magnetic field on 
the observed emission using numerical simulations. In J9J we studied a large number 
of shell collisions with different magnetization levels. In IfTTI we focused on a limited 
number of parameters such as the observer viewing angle, (9, the bulk Lorentz factor of 
the slower shell, T r, and the relative Lorentz factor, A g (1 '/ /F R ) - 1, where T/ is the 
bulk Lorentz factor of the faster shell. In the present work we summarize some of the 
results studied therein. 

We describe the common numerical setup in Sec.0 The results are shown in Sec.0 
and discussed in Sec.0 
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2 Numerical Setup 

A modified version of the code SPEV |J8) was employed to calculate the non-thermal 
emission from ISs. We assume a cylindrical shell geometry and perform all the calcula¬ 
tions in the rest frame of the shocked fluid (RFSF). The shell interaction was simplified 
as a one-dimensional Riemann problem and focus our resources on a more detailed 
treatment of the non-thermal radiation. Our method consists of three steps that we sum 
up in the following paragraphs. 

Solution of the Riemann problem. Employing an exact relativistic magnetohydro¬ 
dynamics (RMHD) Riemann solver we determine the properties of the internal shock 
waves. We follow the procedure described in 13 to set-up the shells and to extract the 
information needed for the steps 2 and 3. 

Non-thermal particles transport and evolution. The particles are injected behind 
the shock fronts following the prescription of Him We assume that a fraction of 
the thermal electrons are accelerated into a power-law distribution at high energies, and 
that their energy density is a fraction of the internal energy density of the shocked fluid. 
In the RFSF the shocks are propagating away from the initial discontinuity, injecting 
and leaving non-thermal particles behind. We evolve the energy distribution of non- 
thermal electrons taking into account synchrotron and inverse-Compton (IC) losses. 
See El for more details. 

Radiative transfer. The total emissivity at each point is assumed to be a combina¬ 
tion of the following emission processes: (1) synchrotron radiation, (2) IC upscattering 
of an internal (SSC) and an external radiation field (EIC). 


3 Results 

We group our models according to the initial shell magnetization, <r := B 2 /4npT 2 c 2 . 
Hereafter the subscripts L and R will denote left (faster) and right (slower) shells, re¬ 
spectively. We compute the spectra for a typical source located at z — 0.5. In the rest of 
this section we will present some of the final SEDs resulting from our simulations. A 
larger collection is shown in G31. The SED of each model has been averaged over the 
time interval 0 - 10 6 s. 

3.1 Weakly-magnetized models 

The SEDs computed for the models with ctr- ctr- 10 -6 ,Tr = 10, f) = 5° and varying 
A g, are shown in the upper left panel of Fig. Q] The spectra show that with increasing 
A g the IC component also increases, up to three orders of magnitude. In order to see 
the effects on each emission process, the synchrotron, SSC and EC components for 
A g = 0.5,2.0 are shown as dashed, dot-dashed and dot-dot-dashed lines, respectively. 
While the three spectral components (synchrotron, SSC and EC) have approximately 
the same order of magnitude for A g = 0.5, for A g = 2.0 the SSC is almost two orders 
of magnitude more luminous than the other two. 

3.2 Moderately-magnetized models 

The SEDs of the family of models ctr = ctr = 10 -2 ,Ag =1.0 and 6 - 5° are presented 
in the upper right panel of Fig. Q] As in the previous section, SSC and EC components 
are shown as dashed, dot-dashed and dot-dot-dashed lines, respectively, but now for 
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Table 1: Physical parameters in the RS shocked region for the family of strongly 
magnetized shells with different A g. The bulk Lorentz factor of both shocked re¬ 
gions is denoted by T, while r, B, Q, y\ and yi denote its compression ratio, co¬ 
moving magnetic field, comoving number of electrons injected per unit volume and 
unit time, and lower and upper cutoffs of the injected electrons in the RS (see Eq. 11 
of MA12). Q n - Qx 10 -11 . t' cr := A//(c|/T|) is the shock crossing time, where A/ 
and /?' are the shell width and the shock velocity in the frame moving with the con¬ 
tact discontinuity separating both shocks (section 2 of MA12). y c := 72/(1 + vo yit' cr ) is 
the cooling Lorentz factor of an electron after a dynamical time scale (shock crossing 
time), vo := ( 4-/3)cctt(u' b + u' ext )/(m e c 2 ) is the cooling term, where <tt is the Thomson 
cross section and the primed quantities are measured in the comoving frame. When 
ydy\ » («) 1 the electrons are slow (fast) cooling. Note that the Q\ 1 for A g = 2.0 is 
much larger than Q\ \ of the other models because y\ - yi. 


the models with Fr = 10,17,25. The synchrotron component for Fr = 10 is — 20 times 
brighter than the SSC one, in contrast to the EC which is 100 times dimmer. For Fr = 25 
the EC is of the same order of magnitude of SSC and synchrotron. The latter two 
decrease one order of magnitude between the models with Tr = 10 and T/,, while the 
EC grows by almost one order of magnitude. This follows from the fact that the number 
of electrons and the comoving magnetic field strength decrease with the increasing T /< 
0, which means that there are less synchrotron photons and less electrons which can 
scatter them in the SSC process. We model the Klein-Nishina cutoff as a sharp cutoff, 
and its effect can be seen at - 10 23 Hz, where all the EC spectra coincide. From the 
inset we can see that there is no significant change in the flux of y-photons, although 
the spectral index softens for increasing Fr. 

3.3 Strongly-magnetized models 

The third family consists of the strongly magnetized models where ay. = 1 and <tr = 0.1. 
The corresponding SEDs of this series of models (lower panel of Fig. [T] are analogous 
to the models of Sec. [3T2l As we can see, for Tr = 10 the synchrotron component is 
- 100 times brighter than the IC. For I r = 25 this difference decreases to one order of 
magnitude. Once again, EIC component rises with rising Tr, to the point in which it 
begins to be comparable to the synchrotron component, developing an IC hump. The 
EIC spectra interestect in a single point due to our treatment of the Klein-Nishina cut¬ 
off. In the inset we can see that the flux of y-ray photons does not change appreciably 
in this family of models. Table |T|lists a number of physical parameters in the reverse 
shock (RS) of the present models. We see that the electrons in the RS of the strongly 
magnetized models are fast-cooling. In fact, for A g =1.5 the injected electron spectrum 
is almost mono-energetic. 
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Figure 1: Upper left panel: Averaged spectra resulting from collisions of weakly mag¬ 
netized shells (ctl = ctr = 10 -6 ). Upper right panel: Same as the upper left panel, 
but for moderately magnetized shells ( ctl = ctr = 10 -2 ) and varying F r. Lower panel: 
Same as upper right panel, but for strongly magnetized shells (cry = 1 ,<tr = 0.1). The 
inset in each panel shows the photon spectral slope T p h as a function of the photon flux 
F p h in the y-ray band (see Sec. 3.2). Colors of the points correspond to the line colors 
in the main plot. For selected models, the synchrotron, SSC and EIC contributions 
(dashed, dot-dashed and dot-dot-dashed lines, respectively) are shown. 
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3.4 y-rays spectral slope 

A linear least-squares algorithm is used to deduce the y-ray spectral slope F p i,. Due 
to the fact that we are not modeling the Klein-Nishina part of the spectrum, we only 
performed the calculations of T p h for those models that do not show a large drop-off 
in the photon flux. In Fig. Owe show T p h as a function of the photon flux for energies 
> 0.2 GeV, where F p h is the photon flux for photon energies >0.1 Ge V [jT). We compare 
our results with sources found in 2LAC |Z|, restricting the comparsion to sources with 
0.4 <z< 0.6. In Fig. [2] we see that weakly and moderately magnetized models overlap 
with the observations, with more weakly than moderately magnetized models falling 
within the observed part of the parameter space. 

3.5 Compton dominance Ac 

Another way to correlate magnetization with observed properties can be found by rep¬ 
resenting the Compton dominance Ac as a function of the ratio of IC-to-synchrotron 
peak frequencies vjc/v syn . Models with intermediate or low magnetization occupy a 
range of Ac roughly compatible with observations 0, while the strongly magnetized 
models tend to have values of Ac that are probably uncompatible with those observed 
in actual sources, unless collisions in blazars happen at much larger Lorentz factors 
than currently thought. 

4 Conclusions 

We have shown that the SEDs of FSRQs and BL Lacs strongly depend on the magne¬ 
tization of the emitting plasma. Our models predict a more complex phenomenology 
than what is currently supported by the observational data. A conservative conclu¬ 
sion is that the observations restrict the probable magnetization of the colliding shells 
that take place in actual sources to, at most, moderate values (i.e., cr < 10 -1 ), and if 
the magnetization is large, with variations in magnetization between colliding shells 
which are smaller than a factor ~ 10. We find that FSRQs have observational prop¬ 
erties corresponding to models with negligible or moderate magnetic fields. BL Lacs 
with moderate peak synchrotron frequencies v syn < 10 16 Hz and Compton dominance 
parameter 0.1 > Ac > 1 display properties that can be reproduced with models with 
moderate and uniform magnetization (cr^ = ctr = 10 -2 ). We find that a fair fraction of 
the blazar sequence can be explained in terms of the intrinsically different magnetiza¬ 
tion of the colliding shells. We observe that the change in the photon spectral index 
(r ph ) in the y-ray band can be a powerful observational proxy for the actual values 
of the magnetization and of the relative Lorentz factor of the colliding shells. Values 
r p h > 2.6 result in models where the flow magnetization is cr ~ 10 -2 , whereas strongly 
magnetized shell collisions (cr > 0.1) as well as weakly magnetized models may yield 

r ph < 2.6. 
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Figure 2: Left panel: Spectral slope Y p h for the photon energies > 200MeV as a func¬ 
tion of the photon flux for energies >100 MeV (T|. The symbols joined by lines 
represent our numerical models, while cyan and magenta triangles represent BL Lacs 
and FSRQs at redshift z — 0.5 from 2LAC 12]. Right panel: Compton dominance Ac 
as a function of the synchrotron peak frequency v syn for all the models studied in HQ. 
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